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Abstract 

Higher-order hadronic contributions to the anomalous magnetic moment of the elec- 
tron, muon, and r lepton are considered in detail. As a main result we find a reduc- 
tion by —11 x 10~ n for the g-2 of the muon as compared to previous calculations. 
Analytical expressions for the kernel functions of higher-order hadronic effects are 
presented. We employ the method of asymptotic expansions to calculate kernel 
functions analytically in terms of a series expansion in the lepton mass. 
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ftp at ttpux2.physik.uni-karlsruhe.de (129. 13. f 02.139) as /ttp96-26/ttp96-26.ps or via www at 
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The investigation of the anomalous magnetic moment of the muon has received con- 
siderable interest recently due to the upcoming BNL experiment [|I]]. The main goal 
of this experiment is the confirmation of electroweak loop effects as predicted by 
the Standard Model. At the same time, deviations may be an indication of physics 
beyond the Standard Model. It is therefore necessary to fix the Standard Model 
prediction as precisely as possible. 

In this context we note that not only electroweak loop effects but also higher-order 
hadronic contributions come into the range of experimental detection. The designated 
experimental uncertainty of ±40 ■ 10 -11 in the Brookhaven experiment allows for an 
establishment of higher-order hadronic effects by at least 4 standard deviations. For 
the electron the situation is similar: from improvements of the Penning-trap tech- 
nique, experimentalists expect a reduction of the present experimental uncertainty 
from Aa e = ±4.0-10~ 12 to Aa e = ±0.6- 10~ 12 |2|, in which case hadronic contributions 
to a e = (g e — 2)/2 would be a 2a effect. 

Present experimental limits for the r lepton are of the order of |Aa r | ~ 6 • 10~ 3 []3], f|], 
nearly three orders of magnitude below the expected hadronic effects. 

The purpose of the present work is to calculate higher-order hadronic diagrams 
for the electron, muon, and r lepton. A new method for the computation of kernel 
functions is presented. Combined with previous work we now provide a complete set 
of analytic kernel functions that allow a numerically reliable evaluation of higher-order 
hadronic contributions to the g-2 of leptons. 

This work is based on a recent analysis of experimental data by Eidelman and 
Jegerlehner We follow this reference in the treatment of integration over hadronic 
data. In particular, we also use the trapezoidal rule for the integration over the 
p resonance. Higher resonances are included using the narrow width approxima- 
tion. The set of kernel functions presented here can, however, also be used for other 
parametrizations of hadronic data. 

The relevant Feynman diagrams are depicted in fig. (2). As in the case of the 
leading order hadronic correction (fig. 1), these diagrams cannot be computed in 
QCD perturbation theory. Instead, experimental data on the branching ratio R(s) 
for e + e~ — ► hadrons can be used to compute the diagrams of figs. (1) and (2a-c). 
This is not possible for the so-called light-by-light-diagram (fig. 3). The latter has 
been computed recently for the muon by two groups using the extended Nambu- Jona- 
Lasinio model || 0]. It will not be discussed in the present paper. 

The paper is organized as follows: In the first section we discuss the calculation 
of kernel functions and collect relevant formulae. In the following sections we apply 
these results to the computation of higher-order hadronic contributions to the muon, 
electron, and tau, respectively. Whenever possible, different methods (analytical and 
numerical) have been used as a cross-check. 
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Kernel Functions 



The common feature of all hadronic contributions discussed in this paper will be 
that they are self-energy insertions in the photon propagator. We introduce the 
polarization function IV for the hadronic insertion with 

IV (g) = % (q^q v - g^q 2 ) I% 2 ) (1) 

ntf) = -I f" ^Jpfl . (2) 

7T J Aral S q 2 — S 

The substitution rule from a photon propagator to a renormalized propagator with 
hadronic insertion follows: 

-^.ir^L^i) . (3) 

q 2 7T J Ami s \ q — s ) 

The calculation of higher-order kernel functions reduces to the computation of 
two-loop QED diagrams where one photon is massless and the other photon has the 
effective mass y/s. The s-integration, which has been shifted to the very end, is 
performed over experimental data using 

a(e + e- -> hadrons) 

R{s) = — — — - = 12vr M s 

cr (e+e — > ) 

and the Born cross section <7o(e + e~ — > /U + /i~) = ^f-- The leading-order hadronic 
contribution (fig. 1) to the anomalous magnetic moment of a lepton a; is given by 



a i 



(i) 



1 fa\ 2 , R(s)KW(s 



3 \7T/ J Ami S 

with 



ds^^ ^ (4) 



where m denotes the mass of the external fermion; explicit expressions for ifW(s) 
can be found e.g. in |5|, |8|]. After renormalization, the higher-order kernel functions 
for the diagrams in fig. (2a,b) can be defined in a similar way: 

' 3 \irj JAml s y ' 

Once the functions K <k2 \s) are known analytically, the higher-order hadronic con- 
tributions to the g-2 of leptons can be computed on an equal footing with the leading 
order (e.g. with the same systematic error). 

The present situation concerning K^ 1 is the following: For the class of diagrams 
in fig. (2a), a complete analytical result for the kernel function has been given in 
(eq.(3.21)). In that reference an expansion to first order in m 2 / s has also been 
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provided. This expansion has been employed for the calculation of hadronic con- 
tributions from diagrams (2a). In the case of diagrams (2b) and (2c) the integral 
representations eq.(§, [13]) have been used. 



As a new approach we employ the method of asymptotic expansions [[L0| . This tech- 
nique has recently been applied to the calculation of two-loop electroweak corrections 
to the g-2 of the muon [[IT], 0. The calculation of hadronic kernel functions for the 
muon and the electron can be treated in exactly the same way. 

The result is an expansion in the small parameter m 2 / s, i.e. the ratio of the muon 
(or electron) mass m and the energy y/s of the final dispersion integral over hadronic 
data. This final integration is numerically well-behaved over the whole region of 
integration. 

For the r lepton such an expansion cannot be performed since the integration 
over hadronic data begins with y/s = 2m n <C m T . However, it is possible to derive 
an analytic formula for the diagram (2b) which describes the contribution for the r 
with sufficient accuracy. 



As a first step, we repeat the calculation of the kernel function for the sum of the 14 
diagrams in fig. (2a) using the asymptotic expansions method and confirm the result 
of [fj. The expansion up to fourth order reads 
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We now discuss diagram (2b): In the following, m denotes the mass of the external 
fermion, rrif the mass of the fermion in the loop. The result of asymptotic expansions 
including fifth order is 
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and can be inserted directly in eq.(||). 
For diagrams (2b) there are other representations (numerical and analytical) which 
may be used to check the result from asymptotic expansions. 
There are two possible starting points for numerical work: 



,(26) 
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with either |13| 
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Finally, we derive an analytical expression for the kernel function, starting from 
eq.fllOD. Inspection of eq.(U) shows that, for light fermions in the loop, the terms 

suppressed with — £ are numerically negligible. We therefore expand (3 in eq . (|l0l) 

m 2 m 

in and retain only unsuppressed terms. After that, the x integration can be 
performed and we obtain 
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and 2+2 and the dilogarithm defined as follows: 
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This expression is exact up to terms of order 0(— £). In particular, we still have the 

full s dependence. Upon expansion of K^ lytical (s) in ^ we recover the unsuppressed 
part of eq.(||). 



For the diagram (2c) with two hadronic self energy insertions on the photon propa- 
gator the representation is |T^] 

a (2c) = — V r dsa H (s) f°° ds'a H (s') C dx -r J (I-.r) 



16vr 5 a74m| ' Jiml ' Jo [x 2 + (1 - x)^][x 2 + (1 - x)-£s\ 

" 3 >" M .ggM _ (13) 



97T 3 J4m2 

An analytical expression for K^ 2c \s, s') is easily found; setting b = and c = ^ 
we obtain for b ^ c: 



, fr +v /-(4-fe)^ 

ir (2c) (s,s') = ~-b-c 
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2 (6 - c) 2 (6 - c) ^- (4 - 6) 6 



, _ , 21 M c 2 (2-4c + c 2 ) ln( c+ ^ (4 - c)c ; 
v — 2 + c) cr ln(c) v ; K c-y/-(4-c)c' 



, >n . + 1 V , (14) 

2 ( b - c) 2 (6 - c) v /- (4 - c) c 

the corresponding formula for 6 = c is 

if( 2c )( s , s ') = --2c + -(-2 + c-4 1n(c)+3cln(c)) (15) 
/ o , „ 2 , c (12-42c+22c 2 -3c 3 ) ln( c+ ^ ( ' 4+c)c ) 

c(— 2 + 4c — CTJ ( V 7 V c- v /(-4+c)c / 



2(" 4 + c ) 2 (-4 + c) J(-4 + c) c 



Muon 

We take the value for the leading-order hadronic contribution to a M from [TjJ (see 



lq , |1T|1 for other recent analyses) 

a { V = 7023. 5(58. 5)(140.9) • 10~ n . (16) 



Here and in the following the first bracket gives the statistical, the second the system- 
atic error. For the higher-order effects we use the combined statistical and systematic 
error of the leading order as an error estimate. 

Higher-order hadronic contributions to the g-2 of the muon were considered for the 
first time in 1976 |T3 and reevaluated in 1984 [T3] with new input from experimental 



data and a new calculation of the light-by- light diagram. 

We start with a discussion of the contribution from fig. (2a), employing the ex- 
pansion eq.(0). If we take only the first line in this formula (this corresponds to eq. 
(3.25) in [|) we obtain 

a^ a) (first order) = -200(4) ■ 10~ n , (17) 



which is very close to the value of —199 • 10 11 in | 18fl . Including the higher-order 
terms shifts the number to 



G (2a) = _2n(5) . io- 11 . (18) 

We checked the integration over the expanded kernel function against the full 
formula of ref. [f|. The deviation is below 0.3% if at least four orders in the expansion 
are included. It should be noted, however, that the numerical evaluation of eq.(3.21) 
in for large values of s (i.e. in a region where the expansion converges rapidly) is 
less reliable than the expansion itself. 



We now turn to the contribution from the diagrams of fig. (2b). We restrict ourselves 
to the case of an electron in the loop. The contribution of the r loop is suppressed 
by . Using the kernel function of eq.(|8|) we obtain in agreement with JTB| 

af ] = 107(2) ■ IO" 11 . (19) 

Again we notice that at least four or five terms in the expansion must be included 
in order to obtain a stable result. The number decreases from 126 • 10 -11 in first order 
to the value of eq. ([l9|) in fifth order (between 4th and 5th order the value changes by 
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0.2%). The part of eq.(BI) suppressed with turns out to be numerically completely 
negligible. 

The result of eq. (|T^) has been checked using different numerical integration meth- 
ods and the analytical representation of eq. (|i~2"D . Whereas the analytical result coin- 
cides with the expansion, care has to be taken in the numerical integrations. 

Numerically, the representation eq.([T0|) is much less appropriate than eq. (|TT|) . The 
latter is a flat function and the double numerical integration over x and y yields the 
result in eq.([L9|). The representation eq.(|T0"D is a peaked function and although for 
fixed s the error is at worst of the order of 1% it adds up systematically in the s 
integration and reaches 10% in the final result. 



To conclude the discussion for the muon, we give the value we obtain for diagram 
(2c) with two hadronic bubbles on the photon propagator. Using eq. (p~3|) we obtain 

ajf) = 2.7(0.1) • 10~ n . (20) 



This number can be reproduced using the analytical expressions eq . ( |14| , [T5D with 
quadruple precision Fortran. 

We find a total shift in the higher-order hadronic contributions of from [|l8|j 

a (2«+26+2c) = _ 90 ( 5 ) . 1Q -11 ( 21 ) 

to 

a (2«+26+2c) = _ 10 i(6) ■ IO" 11 . (22) 

The difference (—11) • 10~ n is somewhat smaller but still of the same order of magni- 
tude as the experimental uncertainty in the BNL experiment, or approximately 25% 
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of the two-loop electroweak contributions. It should be emphasized that the reason 
for the shift doesn't reside predominantly in the input of new experimental data; 
rather, the high precision of the Brookhaven experiment requires the use of the full 
kernel functions. 

Electron 

Higher-order hadronic contributions to the g-2 of the electron have only recently 
been considered for the first time. Some features are special to the electron (as 
opposed to the muon case) and will be discussed in the following. 
For the electron the calculation of diagram (2a) is very similar to the muon case. It 
is sufficient to take the first term in the expansion of eq. ([?]). 

r , 2oW , ml /23, s 7T 2 223 \ , N 

KV a >(s) = -2— e - — In — + . 23 

V ; s V36 m 2 3 54 J K J 

Here, the contribution from the large logarithm dominates in the whole region of 
integration. We obtain the value 

a (2a) = _ 2 25(0.05) ■ 1(T 13 . (24) 

The diagrams of fig. (2b) with an external electron yield a numerically negligible 
contribution. The first observation is that the diagrams with an electron self-energy 
on the photon propagator are already contained in the set of diagrams (2a). Diagrams 
with a muon or a r lepton on the photon propagator are suppressed by additional 
powers of or respectively, and hence can be omitted. This is a special case 

of the general statement that in diagrams of the type (2b) large logarithms occur 
if the mass of the external fermion is bigger than the mass of the fermion in the 

2 

loop (see, for instance, the leading term 5 In— £ in eq.(|8|)); in the opposite case there 
is a mass suppression. The double hadronic bubble diagram (2c) is of the order 
a (2c) ~ 1.2 • 10~ 19 for the same reason. 

For an estimate of the total higher-order hadronic contribution to a e including light- 
by-light scattering we refer to ref. @: 

af a) + af ] = (-2.25 - 0.122) • 1(T 13 = -0.24 • 10~ 12 , (25) 

where the light-by-light value has been estimated by rescaling the corresponding 
number for the muon. 

The higher-order corrections to a^ ad (leading) [|1J 

= 1.8847(0.0165)(0.0375) • 10~ 12 (26) 

reduce the total hadronic contribution by 13% ; the final number is 

a^ ad (total) = 1.645 • 10" 12 . (27) 
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r lepton 

The leading-order hadronic effect is 

4 1} = 338.30 (1.97)(9.12) • 1(T 8 . (2* 



In [|19| an estimate of the higher-order hadronic contributions to a T was given 

^(higher orders) = -1.2 (2) ■ 1(T 7 . (29) 



The number can be reproduced by rescaling the muon results of with the 

2 

factor ^g- and taking the 1991 value for m T : 



f (2a+26+2c) + a li 9 ht-by-li 9 ht\ = ( _ 9Q + 4g) . 1Q -U . /™\ 2 = _ l 2 . 1Q - 7 



It will be shown in the following that, as in the case of the electron, simple rescaling 
of the muon results doesn't give a correct estimate. 

2 

For the r lepton it is no longer possible to use expansions in — over the whole region 
of integration because the integration starts with s = 4m 2 . The diagrams of fig. (2a) 
can still be calculated with the general kernel function of |J. We obtain the value 

a (2a) = _ 3 . 18 (o.08) • 1(T 8 . (30) 

In the case of the diagrams of fig. (2b) we have to sum over the contributions 
of the electron and the muon in the loop. Here, we apply two different methods to 
compute the anomaly: a two-dimensional numerical integration yields the results 

4 26) (electron) = 7.48(0.2) • 10" 8 (31) 
4 26) (niuon) = 1.91(0.05) ■ 10" 8 . (32) 

Using the analytic formula eq.(0) we obtain the following values: 

4 26) (electron) = 7.48(0.2) • 10" 8 (33) 
4 26) (niuon) = 1.88(0.05) ■ HT 8 . (34) 

As expected, the electron contribution is the same and we have a 2% deviation for 
the muon. This is consistent with the order of magnitude of the neglected terms, i.e. 

m l , m>l / \ 

— ^ In ^^ 0.02 (35) 



m 2 . m 2 



For the r lepton, diagram (2c) also yields a numerically important contribution. 
From eq.(|l~3D we get the value 

4 2c) = 1.4(0.1) ■ 10~ 8 . (36) 

The higher-order correction to a^ ad (leading) enhances the total hadronic contri- 
bution by 2% (without the light-by light contribution): 

a (2a+2b+2c) = 7 6 (o.2) ■ 10- 8 . (37) 
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This should be contrasted with the estimate from simple rescaling: 



( 




a 



(2a+2b+2c) 



(-101) ■ 10 



-11 



( 



1777\ 2 



106 J 



2.8 • 10 



,-7 



Keeping this uncertainty in mind, we give an estimate for higher-order hadronic 
contributions to a T with eq.(^) and a rescaled light-by- light value from ||: 



Conclusions 

We presented analytic formulae for the calculation of higher-order hadronic contri- 
butions to the g-2 of leptons and introduced the method of asymptotic expansions as 
a convenient tool to calculate the corresponding kernel functions. We demonstrated 
that in order to match the high experimental precision of the upcoming experiments, 
it is necessary to include the full kernel functions in theoretical predictions. The 
higher-order hadronic effects on the anomalous magnetic moments of the muon, elec- 
tron, and t lepton are presented. The remaining uncertainties in the predictions 
from the Standard Model reside now in the experimental error in the measurement 
of a(e + e~ — > hadrons) and the calculation of light-by-light-scattering effects. 
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Fig. (1): Leading order hadronic contribution 




Fig. (2a): Dashed lines indicate a hadronic insertion on the photon propagator. 
Mirror counterparts and diagrams with interchange of massless and "massive" photon 
propagators have to be included. 




Fig. (2b) Fig. (2c) Fig. (3) 
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